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ABSTRACT 


Analysis  of  data  we  collected  in  south  central  Siberia  and  Northern  Mongolia  along  with  digital 
data  from  the  Global  Seismic  Network  have  been  analysed  to  identify  effects  on  propagation  of  large 
scale  heterogeneity  in  the  mantle  that  need  to  be  taken  into  account  when  using  seismic  arrays  to 
monitor  the  Comprehensive  Test  Ban  Treaty.  We  find  that  in  the  vicinity  of  the  Baikal  rift  zone 
in  Siberia,  P  and  S  waves  exhibit  anomalous  delays  and  attenuation.  S  waves  show  clear  evidence 
of  anisotropy  effects.  This  anomalous  region  that  extends  over  several  hundred  km  either  side  of 
he  rift  zone,  and  is  presumed  to  extend  along  the  600  km  of  its  axis,  represents  an  attenuating 
and  delaying  mantle  barrier  to  seismic  phases.  These  findings  can  be  used  to  infer  properties  of 
analogous  tectonic  regions  world-wide  (Basin  and  Range  of  North  America,  Northern  Tibet,  East 
African  and  Rio  Grande  rifts.  Red  Sea,  Gulf  of  California).  We  have  analysed  all  stations  of  the 
Global  Seismic  Network  (GSN)  for  SKS  splitting  in  the  mantle.  We  compare  SKS  waveforms  with 
normal  mode  synthetics  modified  by  splitting  by  one  or  two  layers  of  anisotropy  in  the  mantle. 
Splitting  is  recognised  as  scattering  of  radial  energy  into  tangential  energy.  We  find  that  stations 
can  be  characterised  by  a  single  layer  of  splitting  in  their  vicinity  but  for  some  arrivals  scattering 
deeper  in  the  mantle  either  due  to  a  second  layer  of  anisotropy  or  3D  scatterers  gives  anomalous 
tangential  energy.  This  observation  has  led  to  development  of  a  new  method  to  identify  deep 
mantle  scatterers,  by  arrivals  with  anomalous  tangential  energy.  We  have  examined  P-S  scattering 
from  12  Lop  Nor  nuclear  explosions  recorded  at  GSN  stations  in  Asia.  We  have  identified  PcS 
and  possible  source  generated  S.  Both  show  splitting  effects  that  must  be  deconvolved  before  they 
can  be  confidently  be  interpreted  as  nuclear  explosion  P-S  scattering.  Lowering  the  threshold  for 
discrimination  and  yield  estimates  of  nuclear  explosions  requires  that  these  effects  on  propagation 
of  the  generated  phases  are  recognized,  and  be  taken  into  account  quantitatively. 

INTR.ODUCTION 


Most  nuclear  tests  occur  on  continents.  Continental  geology  can  broadly  be  separated  into  stable 
shields  and  mobile  belts.  Deformation  of  mobile  belts  gives  rise  to  orogenic  zones  and  rift  zones. 
Recently  it  has  been  recognized  that  large  variations  in  the  seismic  properties  of  the  upper  mantle 
occur  beneath  the  mobile  belts  of  the  continents  and  these  need  to  be  taken  into  account  when 
using  seismograms  to  discriminate  explosions  from  earthquakes,  to  measure  seismic  yield  and  to 
locate  events. 
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This  report  has  three  parts.  Part  one  (sections  1-4)  documents  lateral  variation  in  attenuation, 
anisotropy,  and  seismic  velocity  in  a  1200x400  km  area  centred  on  the  Baikal  rift  zone  in  Siberia,  and 
the  Mongolian  fold  belt,  based  on  two  portable  array  experiments  we  carried  out  with  colleagues 
from  University  of  Wisconsin  and  the  Institute  of  the  Earth’s  Crust,  Irkutsk  in  the  summers  of 
91  and  92.  It  repeats  much  of  what  appears  in  our  previous  DARPA  report  (Davis  and  Gao, 
1996,  DARPA  Final  Report)  but  includes  new  analyses  of  S-wave  delays  and  recognition  of  lower 
mantle  velocity  anomalies  in  the  data  and  possible  slab  effects  in  the  source  zone  of  the  measured 
teleseisms..  Part  two  (section  5)  deals  with  application  of  normal  mode  synthetics  to  scattering  in 
the  mantle  and  SKS  splitting  and  is  a  result  of  a  collaboration  with  Prof.  John  Woodhouse  (who 
has  developed  the  synthetics  code)  while  the  PI  was  on  sabbatical  in  Oxford.  Part  three  is  (section 
6)  reports  on  our  (Laura  Borton,  AASSERT  graduate  student)  attempts  to  recognise  teleseismic  S 
body  waves  from  Lop  Nor  Nuclear  explosions  and  what  effects  the  anomalies  that  we  have  measured 
in  the  previous  sections  have  on  seismograms  from  nuclear  explosions  recorded  by  the  GSN. 

1.  The  Baikal  Rift  Zone  Experiments 

Analysis  of  the  data  from  the  Baikal  rift  zone  experiment  has  revealed  the  distribution  of  attenuation 
and  velocity  anomalies  beneath  the  rift  zone  which  abuts  the  Siberian  craton  to  the  NW  and  the 
Mongolian  fold  belt  to  the  SE.  The  attenuation  analysis  finds  that  relative  and  t*  have  a  0.2  and  5 
second  anomalies  in  the  rift  zone.  Explosions  detonated  in  this  area  would  have  reduced  body  wave 
amplitudes  at  arrays  at  teleseismic  distances  such  as  NORESS  with  a  concomitant  underestimate 
in  yield.  We  have  observed  S  wave  splitting  of  on  average  1  s  at  virtually  every  station,  implying 
that  significant  anisotropy  is  present  in  the  uppermost  mantle.  Though  fast  directions  are  spatially 
coherent  in  local  regions  there  is  significant  variation  region  to  region.  P  wave  travel  time  delays 
show  an  unusual  pattern  of  a  central  peak  of  about  1  sec  surrounded  by  two  troughs  of  about 
0.5  secs  resulting  in  anisotropic  variation  dependent  on  incidence  and  azimuth.  We  attribute  the 
combination  of  low  Q,  anisotropy  and  travel  time  anomalies  to  small  scale  convective  upwelling 
beneath  the  rift  zone.  Comparison  with  similar  experiments  and  the  global  Earth  models  is  taken 
as  evidence  that  these  effects  are  present  world-wide  beneath  extension  zones  on  the  continents 
and  should  be  taken  into  account  as  part  of  any  global  seismic  monitoring  system.  Our  portable 
array  spans  two  stations  which  will  be  used  to  monitor  the  CTBT:  Talya,  Siberia  and  Ulan  Batur, 
Mongolia.  Our  characterisation  of  regional  propagation  effects  in  their  vicinity  will  thus  be  useful 
for  station  calibration. 
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The  data  tapes  from  these  experiments  have  been  submitted  to  the  IRIS  data  center  and  are 
available  to  any  one  who  is  interested  in  using  them.  The  following  sections  introduce  the  Baikal 
rift  1991-92  experiment  and  briefly  describes  the  major  scientiflc  flndings  obtained  using  the  data 
set. 


2.. THE  1991  EXPERIMENT 

The  1991  experiment  took  place  between  July  7  and  October  5.  Twenty-six  seismic  stations  (10 
PASSCAL  Reftek  and  16  U.  Wisconsin  recorders)  were  installed  about  Lake  Baikal  (Figure  2.1 
and  Table  2.1).  The  UCLA  field  trip  team  installed  and  operated  the  10  Reftek  stations.  All 
the  stations  were  equipped  with  1  Hz  3-component  sensors  and  five  of  the  stations  were  co-sited 
with  broadband  (Guralp  and  STS2)  seismometers.  The  seismographs  synchronize  internal  clocks 
to  signals  from  the  Omega  navigation  system  (locked  to  either  Norway  or  Japan)  which  ensured 
that  the  timing  error  for  most  of  the  data  was  less  than  20  ms. 

Station  Information 

•  Recorders-  The  10  Reftek-72A  recorders  were  lent  to  the  experiment  from  the  IRIS/PASSCAL 
Instrument  Center  at  LDEO. 

•  Seismometers-  3  kinds  of  seismometers  were  deployed:  1)  L4C  Mark  Products.  These  are  1  Hz 
3  component  sensors;  2)  Four  Guralp  VBB  CMG3  seismometers  (vertical  components  only)  with 
velocity  response  flat  from  30  sec  to  10  Hz  at  stations  01,  05,  07  and  10;  3)  One  Strekiesen  STS-2 
broad-band  3  component  seismometer. 

For  all  the  10  stations  channel  1  to  3  recorded  the  L4C  seismometers  with  the  vertical  on  channel 
1,  north-south  component  on  channel  2,  and  east-west  component  on  channel  3.  All  the  Guralps 
were  recorded  on  channel  4.  For  the  STS-2  station,  channel  4  was  recorded  vertical  component  of 
STS-2,  channel  5  north-south  and  channel  6  the  east-west. 

•  Station  Locations-  All  the  station  co-ordinates  were  determined  by  using  a  GPS  Pathfinder 
receiver  (borrowed  from  UCLA’s  GPS  group).  The  errors  of  the  measurement  on  horizontal  co¬ 
ordinates  were  less  than  50  meters.  Due  to  its  poor  accuracy  on  determining  elevation,  Russian 
1:200,000  maps  were  used  to  obtain  elevation  instead  of  using  the  GPS  measurement.  Table  2.1 
lists  station  co-ordinates  and  other  information. 

•  DAS  Movement-  Some  of  the  recorders  (i.e.,  DAS,  Reftek  Data  Acquisition  System)  were  moved 
during  the  experiment  for  various  reasons.  The  main  reason  for  the  movement  was  to  replace  a 
malfunctioning  DAS  with  a  better  one  so  as  to  get  some  useful  data  from  each  site.  Detailed 
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information  about  DAS  movement  can  be  found  in  Table  2.2. 

•  Station  Performance  and  Problems-  Installation  and  servicing  were  carried  out  using  a  combination 
land  vehicles,  Aeroflot  helicopters  and  boats  on  lake  Baikal.  Stations  were  serviced  once  every  20 
days  for  3  months.  Extremely  tricky  roads  and  shortage  of  gasoline  limited  the  frequency  6f  services. 
Cloudy  and  rainy  weather  in  the  short  summer  of  Siberia  caused  some  recorder  power  failure.  All 
of  Irhese  conditions,  together  with  the  non-ideal  performance  of  the  Reftek’s  made  a  flnal  rate  of 
success  78%  for  the  10  stations. 


Data  Information 

•  Recording  Parameters-  For  most  of  the  stations  we  programmed  a  10  samples  per  second  (sps) 
continuous  data  stream  to  record  teleseismic  events;  a  50  sps  triggering  stream  to  record  high 
frequency  local  and  regional  events;  and  for  those  stations  co-sited  with  broadband  stations,  another 
continuous  data  stream  with  sample  rate  1  sps  and  32  bits  data  words  was  programmed  to  record 
surface  waves.  It  must  be  pointed  out  that  at  the  beginning  of  the  project,  a  shorter  recording 
length  and  a  higher  trigger  ratio  for  the  triggered  data  stream  were  used  rather  than  the  ’standard 
parameters’  as  shown  on  Table  2.3.  Detailed  parameters  can  be  found  in  the  log  files. 

•  Data  format-  Data  dumped  from  the  station  disks  were  converted  to  SEGY  format  by  using  the 
REF2SEGY  routine  (Version  91.099)  written  by  Early  et  al.  at  the  PASSCAL  Instrument  center. 
The  resultant  SEGY  format  data  can  easily  be  converted  to  other  formats  such  as  SAC,  AH,  and 
SIERRA  by  using  routines  provided  by  the  instrument  center. 

•  File  name  convention  and  data  organisation-  Files  were  sorted  into  sub-directories  according  to  the 
starting  times  and  data  stream  that  they  belong  to.  For  instance,  under  sub-directory  R188.02  are 
files  in  data  stream  02  with  starting  time  within  Julian  day  188. 

File  names  contain  the  information  about  starting  hour,  minute,  and  second  of  the  file,  as  well 
as  DAS  serial  numbers  and  channel  numbers.  For  instance,  file  20.45.55.0369.2  under  sub-directory 
R188.01  starts  between  188:20:45:55.00  and  188:20:45:56.00  (the  exact  starting  time  of  each  file  can 
be  found  in  the  header  of  the  file);  the  file  was  recorded  by  DAS  369  which  was  station  01  (Table 
2.2)  on  channel  2  (i.e.  the  north-south  component). 

Earthquake  Information 

•  Event  Statistics-  According  to  the  NOAA  Earthquake  Bulletin,  314  events  with  mb  >  5.0 
occurred  during  July  7  and  October  5,  1991  (Figure  2.2  and  Table  2.4).  54  of  them  are  larger  than 
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164  of  these  314  events  range  from  30°  to  85°  to  the  center 
with  mb  >  4.0  occurred  in  the  area  with  delta  distance  less 
events  with  depth  larger  than  100  km. 

Example  of  Seismograms 

''Figures  2.3  and  2.4  show  some  seismograms  recorded  by  the  profile.  Figure  2.3  compares 
teleseisms  and  spectra  from  a  Northern  California  event  fi-om  the  4  Guralp  stations.  Figure  2.4  is 
a  seismic  section  from  a  teleseismic  event  from  1  Hz  vertical  sensors,  on  data  stream  1,  i.e.  the 
continuous  stream. 


mb=5.5.  The  epicentral  distances  for 
of  the  profile.  As  many  as  111  events 
than  30°.  Also  recorded  were  82  deep 


Participating  Personnel 

Personnel  involved  in  the  installation  and  operation  of  the  10  stations  include: 

•  UCLA:  P.M.  Davis,  C.  Davis,  S.  Gao,  H.Liu,  and  P.D.  Slack 

•  Institute  of  Earth’s  Crust,  Irkutsk:  Yu.  A.  Zorin,  A.  Masalski,  V.M.  Kozhevnikov,  V.V.  Mordvi- 
nova,  and  T.  Perepelova.  Academician  N.  A.Logatchev,  director  of  the  Institute,  performed  excellent 
leadership  to  the  entire  project. 

•  Institute  of  Physics  of  the  earth,  Moscow:  M.  Kogan,  S.  Panasyuk,  and  M.  Karpechov. 

•  PASSCAL  Instrument  Center,  LOGO:  R.W.Busby 

In  addition,  our  colleagues  from  U.of  Wisconsin,  Madison  led  by  Prof.  R.P.  Meyer  helped  in 
the  operation  of  the  10  stations  that  were  primary  responsibility  of  the  UCLA  team. 

3.  THE  1992  EXPERIMENT 

During  June  11  and  October  5, 1992  the  second  phase  of  the  experiment  was  conducted.  The  UCLA 
field  trip  team  installed  and  operated  the  twenty-eight  Reftek  stations  and  two  UW  recorder  stations 
(Figure  3.1).  All  the  stations  were  equipped  with  3-component  short  period  sensors  with  central 
frequencies  range  from  0.5  to  2  Hz  and  eleven  of  the  stations  were  co-sited  with  broadband  (Guralp 
and  STS2)  seismometers. 


Station  Information 

•  Recorders-  The  Reftek-72A  recorders  were  lent  to  the  experiment  from  the  IRIS/PASSCAL 
Equipment  Center  at  LDEO. 
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•  Seismometers-  Six  kinds  of  seismometers  were  deployed:  1)  Thirteen  L4C  1  Hz  three  separate 
component  sensors.  2)  Eleven  S13  IHZ  three  separate  component  sensors  from  the  IRIS/PASSCAL. 
3)  Four  L22  2HZ  three  component  sensors  from  the  IRIS/PASSCAL.  4)  Four  Guralp  VBB_CMG3s 
(vertical  components  only)  with  velocity  response  flat  from  30  sec  to  10  Hz  at  stations  06,  15,  80, 
84  and  87  (Table  3.1  and  Figure  3.1).  5)  Five  Strekiesen  STS-2  broad-band  3-component  sensors 
co-*bccupied  with  short  period  sensors  at  stations  00,  01,  11,  16,  81  and  90.  6)  Two  UW  recorders 
along  with  two  HSIO  IHz  sensors  from  UW  were  installed  at  station  08  and  10. 

For  all  the  stations  with  short  period  sensors  the  vertical  component  was  on  channel  1,  north- 
south  component  on  channel  2,  and  east-west  component  on  channel3.  All  the  Guralps  were 
recorded  on  channel  4.  For  the  STS-2  stations,  channel  4  was  recorded  vertical  component,  channel 
5  north-south  and  channel  6  east- west. 

•  Station  Locations-  The  same  method  as  what  was  used  in  the  1991  experiment  was  used  to 
determine  station  locations. 

•  DAS  Movement-  Detailed  information  about  DAS  movement  can  be  found  in  Table  3.2. 

•  Station  Performance  and  Problems-  Stations  were  serviced  once  about  every  20  days  for  about  4 
months.  The  rate  of  success  is  89%  for  the  30  stations. 

Data  Information 

•  Recording  Parameters-  For  most  of  the  stations  we  programmed  a  10  sps  continuous  data  stream 
to  record  teleseismic  events;  a  25  sps  (for  STS2  stations)  or  50  sps  triggering  stream  to  record 
high  frequency  local  and  regional  events;  and  for  those  stations  co-sited  with  broadband  stations, 
another  continuous  data  stream  with  sample  rate  1  sps  and  32  bit  data  words  was  programmed 
to  record  surface  waves  with  wide  dynamic  range.  See  Table  3.3  for  the  ’’standard  parameters”. 
Detailed  parameters  can  be  found  in  the  log  flies. 

•  Data  format-  The  same  as  in  the  1991  experiment. 

•  File  name  convention  and  data  organisation-  The  same  as  in  the  1991  experiment. 

Earthquake  Information 

•  Event  Statistics-  According  to  the  NOAA  Earthquake  Bulletin,  464  events  with  mb  >  5.0 
occurred  during  June  11  and  October  5,  1992  (Figure  3.2).  113  of  them  are  larger  than  mb=5.5. 
The  epicentral  distances  for  217  of  these  464  events  ranges  from  30°  to  85°  to  the  center  of  the 
profile.  As  many  as  191  events  with  mb  >  4.0  occurred  in  the  area  with  delta  distance  less  than 
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30°.  Also  recorded  were  81  deep  events  with  depth  larger  than  100  km. 


Example  of  Seismograms 
Figures  3.3  shows  a  typical  record  section. 

•s. 

Participating  Personnel 

Personnel  involved  in  the  installation  and  operation  of  the  30  stations  include: 

•  UCLA:  P.M.  Davis,  S.  Gao,  H.  Liu,  RD.  Slack,  M.  Benthien,  and  Don  Daniels 

•  Institute  of  Earth’s  Crust,  Irkutsk:  Yu.A.  Zorin,  A.  Masalski,  V.M.  Kozhevnikov,  V.V.  Mordvinova, 
and  T.  Perepelova. 

•  PASSCAL  Instrument  Center,  LOGO:  R.W.  Busby. 

4.  DATA  ANALYSIS  AND  RESULTS 

Using  the  unique  data  set  we  investigated  various  properties  of  the  upper  mantle  beneath  the 
Baikal  rift  zone.  These  include  1).  travel  time  residuals;  2).  seismic  anisotropy;  and  3)  seismic 
anelasticity.  Detailed  descriptions  of  these  studies  are  in  several  published  papers  and  manuscripts 
in  preparation.  In  this  section  we  only  briefly  introduce  some  major  results. 

P-wave  travel  time  residuals  for  the  1991  array 

Upper  mantle  velocity  structure  is  crucial  for  phase  identification  and  location  for  regional  and 
teleseismic  events.  A  three  dimensional  downward  projection  inversion  method  is  used  to  invert  the 
P-wave  velocity  structure  with  constraints  from  deep  seismic  sounding  data.  Our  results  using  the 
1991  data  set  suggest  that  [Gao  et  al.,  1994a]  (1)  the  lithosphere-asthenosphere  transition  upwarps 
beneath  the  rift  zone.  (2)  the  upwarp  has  an  asymmetric  shape;  (3)  the  velocity  contrast  is  -4.9% 
in  the  asthenosphere;  (4)  the  density  contrast  is  -0.57%. 

For  the  inversion  we  used  638  travel  times  from  66  teleseismic  events  at  22  stations  (10  UCLA 
and  12  UW)  and  a  polynomial  up  to  the  6th  order  to  represent  the  lithosphere-asthenosphere 
transition'  The  square  root  of  the  unmodeled  variance  is  0.11  seconds,  which  is  larger  than  the 
observational  error  range  but  is  presumed  to  arise  from  over-simplification  of  the  structure  as  a 
single  undulating  boundary  separating  two  regions  of  constant  but  different  velocity  (Figure  4.1a). 
The  uncertainty  in  the  depth  of  the  boundary  corresponding  to  the  travel  time  misfit  is  ±17  km. 
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The  results  of  the  inversion  are  shown  in  Figures  4.1a  and  4.1b.  The  depth  of  the  interface  ranges 
from  93  km  at  300  km  west  of  the  lake,  beneath  the  craton,  34  km  beneath  the  lake,  and  back  to 
49  km  at  200  km  east  of  the  lake,  beneath  the  fold  belt.  The  upwarp  is  asyminetrical  in  terms  of 
the  difference  between  the  slopes  of  the  western  and  eastern  flanks.  The  western  edge  is  steeper. 
The  asymmetrical  shape  is  consistent  with  the  results  of  gravity  (Figure  4.1b)  and  magnetotelluric 
stiMies  [Zorin,  1971;  Popov,  1990].  The  highest  point  of  the  lithosphere-asthenosphere  boimdary  as 
constrained  by  DSS  data,  reaches  34  km  hear  the  east  shore  of  the  lake,  in  a  region  of  hot  springs. 

P-  and  S-wave  travel  time  residuals  from  1992  long  line  array 

The  data  from  the  array  of  Figure  3.1  consist  of  2128  P-wave  travel  times  from  155  teleseismic 
events,  corrected  by  subtracting  theoretical  arrival  times  [from  the  lASPEI  1991  Earth  model, 
Kennett  and  Engdahl,  1991].  Relative  residuals  were  formed  for  each  event  by  subtracting  the 
event’s  mean  residual  from  the  raw  residuals  and  detrended  to  isolate  velocity  anomalies  having 
wavelengths  less  than  the  array  length  (1260  km). 

To  examine  the  event  location  dependence  of  the  relative  travel  time  residuals,  we  group  the 
events  into  26  groups  by  azimuth  (</>)  and  epicentral  distance  (A)  of  the  sources.  This  gave  517 
mean  travel  time  residuals  from  26  groups  which  were  then  used  in  an  inversion  for  structme. 

Most  of  the  event  clusters  recorded  by  the  1992  profile  (Figure  3.1  and  4.2a)  display  a  peak 
in  their  travel  time  residual  curves  approximately  in  the  region  -30  km  <  a:  <  60  km  with  valleys 
on  either  side.  The  peaks  are  interpreted  to  arise  from  upwarped  low-velocity  structure  associated 
with  convective  upwelling  in  the  mantle.  The  relatively  uniform  distance  between  the  valleys,  and 
the  large  shift  in  their  location  as  a  function  of  incoming  azimuth  axe  interpreted  to  be  caused  by 
two  high  velocity  structures  located  about  ±200  km  deep. 

We  used  an  isotropic  mantle  model  to  interpret  the  travel  time  data  [Gao,  1995].  The  model  has 
a  curved  lithosphere/asthenosphere  boundary  that  upwarps  beneath  the  rift  and  downwarps  on  the 
flanks.  At  the  center  the  asthenosphere  upwarps  from  200  km  to  about  45  km.  The  magnitudes  of 
the  downwarps  are  about  120  km  in  the  west  and  40  km  in  the  east.  The  velocity  contrast  between 
the  lithosphere  and  the  asthenosphere  is  about  -2.4%  (Figure  4.2b). 

Figure  4.2c  shows  mean  S-wave  relative  travel  time  residuals  from  4  events  from  the  South.  A 
peak-to-peak  anomaly  of  about  4  s  is  observed.  The  location  of  the  maximum  anomaly  shifts  about 
70  km  towards  the  NW  from  the  center  of  Lake  Baikal. 
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Seismic  anisotropy  and  mantle  flow 

Flow  in  the  mantle  is  thought  to  orient  olivine  crystals,  resulting  in  anisotropy  in  seismic 
velocities  [Babuska  and  Cara,  1991;  Silver  and  Chan,  1991;  Vinnik  et  al,  1992;  Makeyeva  et  al, 
1992]  with  the  fast  direction  oriented  in  the  direction  of  flow.  Horizpntal  anisotropy'  has  been 
measured  on  the  axis  of  the  Rio  Grande  rift  [Sandvol  et  al,  1992]  but  at  present  measurements 
ha^  not  been  made  on  its  flanks.  Twenty-eight  digital  seismic  stations  and  35  analog  stations  were 
used.  We  find  that,  like  the  mid-ocean  ridges,  away  from  the  rift  axis  the  fast  direction  is  horizontal 
and  normal  to  the  rift  axis  (Figure  4.3).  The  final  results  show  that  [Gao  et  al.,  1994b;  Gao,  1995]: 

1)  Beneath  the  rift  zone  and  Siberian  platform,  most  of  the  fast  directions  are  northwest- 
southeast,  i.e.  orthogonal  to  the  strike  of  the  rift,  and  parallel  to  the  extensional  stress  direction 
[Zohack,  1992]. 

2)  At  the  transition  to  the  fold  belt  in  northern  Mongolia,  the  fast  direction  changes  to  nearly 
east- west,  i.e.  parallel  to  the  faulting  and  fold  axis. 

3)  The  transition  takes  place  between  stations  20  and  22,  over  a  distance  of  100  km. 

4)  The  splitting  ranges  from  0.3  to  2.1  seconds  which  is  consistent  with  a  layer  of  30  to  170  km 
thick  characterised  by  4%  anisotropy. 

5)  Fast  directions  in  the  inner  region  of  the  Baikal  rift  zone  are  distributed  in  two  orthogonal 
directions,  NE  and  NW,  approximately  parallel  and  perpendicular  to  the  NE  strike  of  the  rift. 

As  suggested  in  Gao  et  al.  [1994b;  1997],  the  observed  rift-orthogonal  orientations  of  the 
anisotropy  could  be  the  result  of  uppermost  mantle  flow  associated  with  recent  tectonics.  The 
southern  results  can  be  explained  by  the  shear  field  associated  with  the  collision  of  India  and 
Asia.  The  northern  results  parallel  the  opening  of  lake  Baikal.  From  the  travel  time  anomalies 
we  have  inferred  that  the  asthenosphere  upwarps  in  a  broad  region  extending  either  side  of  Lake 
Baikal,  reaching  the  base  of  the  crust  at  the  lake  itself.  The  inferred  flow  responsible  for  the 
upwarp  is  expected  to  anneal  out  ancient  anisotropy.  Like  the  oceanic  case,  away  from  the  rift, 
the  fast  direction  is  normal  to  the  rift  axis.  The  rift-parallel  fast  directions  near  the  rift  axes  can 
be  interpreted  by  oriented  magmatic  cracks  in  the  mantle  or  small  scale  mantle  convection  with 
rift-parallel  flow.  The  agreement  between  stress  estimates  and  corresponding  crack  orientations 
lends  some  weight  to  the  suggestion  that  the  rift-parallel  fast  directions  are  caused  by  oriented 
magmatic  cracks. 
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Seismic  P-  and  S-wave  attenuation 

We  applied  a  common  spectrum  method  [Halderman  and  Davis,  1991]  to  P  wave  and  S  wave 
data  from  Baikal  to  obtain  relative  t*  values.  The  stations  analysed  were  located  along  the  1280 
km  long  profile  traversing  the  Siberian  platform,  the  BRZ,  and  the  Mongolian  fold  belt’;  and  in  a 
300  by  300  km  cluster  around  the  southern  part  of  Lake  Baikal. 

WVe  used  37  events  for  P-wave  attenuation  studies,  and  5  for  S-wave  studies.  For  the  calculation 
of  both  St*  and  St^,  a  time  window  of  12.8  s  was  taken  starting  approximately  6  s  before  the  onset 
of  the  signal,  and  was  tapered  and  Fourier  transformed.  A  noise  window  of  the  same  length  was 
also  selected  before  the  signal  window  and  Fourier  transformed.  The  noise  spectrum  was  subtracted 
from  the  signal  spectrum. 

The  number  of  measurements  for  St*  and  St^  is  596  and  82,  respectively.  The  mean  values  at 
each  station  and  their  standard  deviations  are  shown  in  Figure  4.4.  The  solid  curves  in  Figure  4.4 
are  obtained  by  averaging  the  measurements  in  a  window  of  200  km  in  width  and  moving  at  10  km 
each  step.  The  peak  to  peak  anomaly  of  0.2  s  for  St^,  and  0.8  s  for  St^  is  observed  in  the  smoothed 
mean  curves.  The  mean  and  standard  deviation  for  Stp  is  0.06±  0.03  s,  and  for  Stg  is  0.14±0.08 
s.  The  Stt  and  Stl  anomalies  associated  with  the  BRZ  are  about  6  times  the  standard  deviation. 
For  St*,  the  peak  values  are  located  in  the  area  of  ±200  km  about  Lake  Baikal,  for  St*,  the  peak  is 
located  at  -160  km  (i.e.,  to  the  Northwest).  Because  all  the  events  for  St*g  were  from  the  south  and 
with  an  S-wave  angle  of  incidence  of  about  35°,  the  northward  shift  may  indicate  that  the  depth 
of  the  center  of  the  low-Q  body  is  located  at  about  240  km. 

Given  that  most  of  the  lateral  variation  of  shear  wave  velocities  are  found  in  the  top  300  km  of 
the  Earth  (e.g.  Woodhouse  and  Dziewonski,  1984),  we  assume  that  the  upper-most  300  km  of  the 
Earth  beneath  the  array  is  responsible  for  the  observed  St*  and  St*.  Under  the  assumption  that 
the  mean  Qp  is  400  beneath  the  Siberian  platform,  the  0.2  s  in  St*  results  in  a  Qp  of  130  beneath 
the  BRZ,  and  the  0.8  s  in  St*g  yields  a  Qs  of  60  beneath  the  BRZ  if  a  Qs  of  200  is  assumed  beneath 
the  Siberian  platform. 

Figure  4.5  shows  the  relation  between  the  smoothed  St*  and  St*.  The  best  fitting  relation  is 
Stl  =  (3.15  ±  0.18)Jt*  —  (0.11  ±  0.01).  The  cross  correlation  coefficient  between  the  two  is  0.83. 
This  ratio  indicates  that  the  observed  St*  is  the  result  of  difference  in  purely  intrinsic  attenuation 
rather  than  difference  in  scattering  (Walsh,  1969;  Solomon  and  Toksoz,  1970;  Richards  and  Menke, 
1983). 


12 


Possible  lower  mantle  velocity  anomalies  and  slab  effects 

High  quality  short-period  data  recorded  by  the  1280  km  long,  50-1-  station  Baikal  network  from 
4  deep  earthquakes  located  in  the  south-west  Pacific  area  (Figures  4.6)  are  used  to  identify  lower 
mantle  velocity  anomalies  and  possible  slab  effects  [Gao  et  ah,  1996].  Th?  optimal  differential  travel 
times  and  t*  between  ScS  and  S,  PcP  and  P,  and  ScP  and  P  are  obtained  by  cross-correlating  one 
ph^e  with  the  waveform  of  the  other  one.  In  total  183  differential  travel  time  residuals  were  used, 
including  96  ScS-S,  47  PcP-P,  and  40  ScP-P  residuals.  Measurements  from  Event  1  are  shown  in 
Figure  4.7.  The  combined  PcP-P,  ScP-P,  and  ScS-S  differential  travel  time  residuals  are  shown  in 
Figure  4.8.  All  the  differential  travel  times  have  a  positive  slope  toward  the  source  direction. 

Some  of  the  possible  explanations  for  the  observations  include: 

1) .  Source  mislocation  The  largest  mismatch  between  the  ISC  and  CMT  catalogs  is  about  15 
km  for  the  events  used.  The  depths  from  the  CMT  catalog  are  always  larger  than  those  from  the 
ISC  catalog  and  this  results  in  a  positive  slope  on  the  travel  time  residuals  toward  the  source. 
However,  the  slope  caused  by  this  effect  is  small.  For  instance,  for  an  event  with  a  depth  of  515 
km,  a  mislocation  in  depth  of  15  km  would  result  in  an  error  on  ScS-S  residuals  of  about  0.9  s  at 
A  =  40°;  and  at  A  =  55°,  it  would  be  about  0.5  s.  Therefore,  the  contribution  to  the  ScS-S  slope 
in  Figure  4.8  would  be  about  -0.053  s/deg,  which  is  about  1/8  of  the  total  observed  slope. 

2) .  Deep  slab  effects  such  as  multipathing  [Silver  and  Chan,  1986]  and  slab  diffraction  [Vidale, 
1987].  The  amount  of  waveform  distortion  and  travel  time  anomalies  depends  on  many  currently 
unknown  facts,  such  as  the  depth  extension  of  the  slabs,  their  dip  angle  and  geometry,  velocity 
distribution  inside  and  outside  the  slabs,  and  the  location  of  the  earthquakes  relative  to  the  slabs. 
At  this  stage  the  influence  of  these  effects  on  the  measured  residuals  cannot  be  quantitatively 
estimated. 

3) .  Lower  mantle  velocity  anomalies.  A  previous  study  by  Wysession  et  aJ.  [1994]  indicates 
that  beneath  South  China  the  D”  shear  wave  velocity  is  high  to  the  north  and  low  to  the  south.  At 
the  north  end  of  the  array  of  bounce  points  (Figure  4.8)  ScS-S  residuals  are  about  -4  s  while  those 
of  ScP-P  are  about  -0.3  s,  implying  that  the  high  shear  wave  velocity  anomaly  is  located  in  the  cS 
branch  of  the  ScS  ray  path.  The  large  ScS-S  negative  residuals  (-4  s)  and  small  PcP-P  negative 
residuals  (-0.3  s)  may  imply  a  decrease  of  the  Poisson’s  ratio  in  the  lower  mantle  high  shear  wave 
velocity  zone. 
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5.  GLOBAL  SKS  ANISOTROPY  AND  SKS  SCATTERING 


Paul  Davis 

Currently  used  seismic  discriminants  and  yield  estimates  rely  on  teleseismic  P  wave  g,mplitude, 
Lg  and  P /Lg  ratios.  Teleseismic  P  waves  are  affected  by  scattering  and  attenuation  in  the  crust  and 
maptle.  A  nuclear  explosion  generates  primarily  a  teleseismic  P  wave.  However  multiple  scattering 
from  interfaces  generate  S  waves  as  well  as  P  waves.  In  order  to  understand  scattering  and  identify 
P-S  scattered  energy  from  major  interfaces  in  the  mantle  it  is  necessary  to  evaluate  the  role  of 
anisotropy  in  the  mantle  in  the  scattering  process.  One  effect  of  anisotropy  on  S  waves  such  as  SKS 
involves  scattering  of  radially  polarised  energy  into  the  tangential  component  as  well  as  changing 
the  relative  phase.  For  SKS  splitting  by  a  single  layer  the  tangential  amplitudes  are  proportional 
to  the  derivative  of  the  radial  amplitudes.  However  for  more  complex  layering  the  relationship  is 
more  involved.  We  examined  SKS  splitting  at  about  150  stations  of  the  global  network  for  the 
years  1998-1995. 

Data  Analysis 

A  filter  was  written  to  extract  all  long  period  SKS  SKKS  seismograms  from  the  Oxford  juke  box 
data  base  at  all  GEOSCOPE  and  GSN  stations  for  the  period  Jan  1988  -  May  1995  corresponding 
to  the  following:  magnitude  greater  than  6;  epicentral  distance  greater  than  90  degrees,  for  a  time 
interval  SKS/SKKS-20  to  SKS/SKKS-l-40  seconds.  This  resulted  in  approximately  100  seismograms 
per  station  per  year.  Normal  mode  synthetics  were  computed  for  each  seismogram  using  Harvard 
moment  tensors  for  source  parameters.  The  normal  mode  program  (Woodhouse  1988)  computes 
seismograms  for  periods  up  to  6  seconds.  Seismograms  and  synthetics  were  rotated  into  radial  and 
tangential  components,  and  in  order  to  extract  good  signal  to  noise  ratio  events,  only  those  for 
which  the  cross-correlation  between  the  observed  (phase  shifted  if  necessary)  and  synthetic  radial 
components  was  greater  than  0.75  were  retained  for  further  analysis.  This  resulted  in  reducing  the 
number  of  useful  seismograms  by  a  factor  of  about  10  for  a  given  station,  but  was  necessary  because 
SKS  arrivals  are  weak  and  often  contaminated  by  microseismic  noise  of  similar  frequencies. 

Standard  1-layer  SKS  splitting  was  performed  on  the  resulting  seismograms.  Seismograms  were 
rotated  into  fast  and  slow  directions  and  the  slow  component  advanced  in  phase,  before  rotating 
back  into  radial  and  tangential  directions.  A  grid  search  of  candidate  fast  and  slow  directions 
was  made  to  find  the  pair  that  gave  minimal  tangential  energy.  Because  the  long  period  data 
is  sampled  at  1  second  intervals  we  perform  phase  shifting  in  the  frequency  domain  which  allows 
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continuous  phase  shifting  compared  to  discrete  shifts  at  the  digitization  interval  in  the  time  domain. 
After  identifying  optimal  parameters  from  the  grid  search  a  non-linear  least  squares  inversion  was 
carried  out  using  the  identified  parameters  as  starting  values.  Final  values  had  formal  errors  of 
several  percent  for  well-determined  minima.  Splitting  parameters  were  determined  for  each  event  at 
station,  but  variability  in  the  results  and  the  means,  depending  on  how  the  results  from  individual 
stations  were  stacked  prompted  us  to  make  a  further  detailed  analysis. 

We  have,  so  far,  examined  in  detail  seismograms  from  GSN  stations  PAS  and  ANMO.  It  appears 
that  the  variability  in  splitting  parameters  could  be  attributed  to  several  causes:  (1)  incident 
polarisation  lying  along  or  within  5  degrees  of  the  fast  direction  or  the  direction  orthogonal  to  it  (the 
null-directions)  gives  rise  to  insuflicient  tangential  energy  to  be  observable  above  the  background 
noise  (2)  scattering  effects  or  anisotropy  effects  not  taken  into  account  in  a  single  layer  analysis 
caused  unmodeled  fluctuations  and  (3)  noise  of  unknown  origin.  We  thus  restricted  our  examination 
to  seismograms  for  which  the  tangential/radial  energy  ratio  was  high,  (greater  than  0.3  to  avoid  the 
null-direction  problem)  and  for  which  the  cross  correlation  between  the  fast  and  slow  components 
was  high  (greater  than  0.8  to  avoid  those  with  significant  scattered  energy).  We  found  that  most 
such  seismograms  could  be  explained  by  a  single  layer  of  azimuthal  anisotropy  lying  beneath  the 
station.  Of  the  64  seismograms  from  ANMO,  19  are  fit  well  by  a  single  layer. 

Having  thus  determined  the  splitting  parameters  for  this  layer  (50  degrees,  1.6  sec.  which 
compares  well  with  previous  estimates  of  40  degrees  and  1.5  sec.)  we  searched  for  effects  of  a 
second  layer.  We  modified  the  normal  mode  synthetics  by  passing  them  through  a  single  layer 
or  double  layer  of  anisotropy  and  compared  them  with  the  radial  and  azimuthal  seismograms. 
Figure  5.1  shows  the  comparison  of  observed  radial  (dotted  curve)  and  synthetic  (smooth  curve) 
seismograms  on  the  left,  azimuthal  components  for  the  one  layer  model  on  the  right  and  those  for 
a  candidate  two  layer  model  in  the  centre.  Event  number,  azimuth,  first  and  second  layer  splitting 
parameters  and  CMT  name  are  listed  on  the  extreme  right  of  the  figure.  Events  23,24,29  and  30 
are  well-described  by  a  single  layer  model.  Event  27  is  better  described  by  the  two  layer  model. 

Figure  5.2  shows  those  events  at  ANMO  that  are  better  described  by  2  layers  of  anisotropy 
and  for  which  the  1  layer  model  gives  a  significant  misfit.  Because  the  second  layer  parameters 
differ  greatly  from  event  to  event  and  because  19  events  are  well  explained  by  a  single  layer,  some 
of  which  have  similar  incident  paths  to  these  events,  we  conclude  that  the  second  layers  can  not 
lie  beneath  ANMO.  Otherwise  all  incident  splitting  would  require  inclusion  of  this  layer  and  one 
second  layer  would  explain  all  the  data.  Anisotropy  deep  in  the  mantle,  perhaps  in  the  D”  region 
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could  explain  these  apparent  second  layer  effects. 

We  are  not  yet  certain  that  the  apparent  second  layer  effects  are  from  another  layer  of  splitting. 
The  fundamental  observation  is  that  a  single  layer  does  not  suffice  and  certainly  a  second  layer 
does  a  better  job,  but  is  this  due  to  more  parameterization  that  allows  a  phase  shift  and ‘amplitude 
change,  but  maybe  is  just  fitting  scattered  energy.  In  either  case  this  is  an  important  effect  which 
is  ih  the  data  at  every  station  and  must  be  taken  into  account  in  comparing  observed  S  waves  with 
synthetics.  Extension  of  this  analysis  to  PAS  reveals  similar  behaviour:  a  single  layer  explains  most 
clearly  well-split  arrivals,  but  a  significant  number  require  a  second  layer  or  significant  modification 
of  the  single  layer  results. 

6.  S  PHASES  FROM  LOP  NOR  NUCLEAR  EXPLOSIONS 


Introduction 

Nuclear  explosions  differ  from  earthquakes  in  that  they  typically  do  not  produce  shear  waves, 
thereby  providing  a  discriminant.  However,  a  series  of  nuclear  explosions  detonated  in  Lop  Nor, 
China  reveal  possible  S  waves.  These  waves  could  be  the  result  of  heterogeneity  in  the  source 
region,  tectonic  strain  release,  or  cavity  collapse.  They  could  also  be  P-to-S  conversions  at  the 
Moho  beneath  the  test  site. 

Data 

The  times,  locations,  and  magnitudes  of  many  of  the  nuclear  explosions  were  obtained  from 
Levshin  and  Ritzwoller  (1995).  Events  occurring  more  recently  came  from  IRIS  (Incorporated 
Research  Institutions  for  Seismology)  by  specifying  parameters  typical  of  explosions  set  off  around 
Lop  Nor:  a  latitude  range  of  40-43  degrees,  a  longitude  range  of  87-92  degrees,  and  a  depth  of 
0-5  km  from  January  1,  1970  through  September  5,  1996.  Table  1  lists  the  twelve  explosions. 
Seismograms  from  stations  with  an  epicentral  range  of  20-100  degrees  were  studied  to  determine  if 
the  Lop  Nor  explosions  generated  any  shear  waves.  No  shear  waves  are  produced  at  distances  larger 
than  100  degrees,  and  crustal  structure  interferes  with  body  waves  less  than  20  degrees.  The  signals 
were  recorded  by  the  Global  Seismic  Network  on  three-component  (North-  South,  East- West,  and 
Vertical)  broadband  seismometers  and  retrieved  from  IRIS.  The  data  was  not  archived  at  all  the 
stations  nor  on  all  the  components;  thus  each  explosion  contains  a  different  number  of  records. 

The  seismograms  in  the  vicinity  of  the  expected  S-wave  arrival  time  were  viewed  using  SAC 
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n 

Latitude 

Longitude 

Depth 

Date 

Time 

Mb 

Ms 

1 

41.75 

88.47 

0 

09/29/88 

07:00:03.1 

4.3 

4.3 

2 

41.57 

88.69 

0 

05/26/90 

07:59:57.8 

3 

41.56 

88.77 

0 

08/16/90 

04:59:57.6 

6.2 

4 

41.60 

88.81 

0 

05/21/92 

04:59:57.5 

.  6.5 

"5.0 

5 

41.76 

88.39 

0 

09/25/92 

0-7:59:59.9 

5.0 

6 

41.67 

88.69 

0 

10/05/93 

01:59:56.6 

5.9 

4.7 

-  7 

41.57 

88.70 

0 

06/10/94 

06:25:58.0 

5.7 

8 

41.55 

89.07 

0 

10/07/94 

03:25:58.0 

5.9 

9 

41.60 

88.82 

0 

05/15/95 

04:05:57.8 

6.1 

10 

41.56 

88.80 

0 

08/17/95 

00:59:57.9 

6.0 

11 

41.66 

88.69 

0 

06/08/96 

02:55:57.9 

5.9 

12 

41.82 

88.42 

0 

07/29/96 

01:48:57.8 

4.9 

Table  1 


(Seismic  Analysis  Code).  For  each  event,  the  three  components  were  plotted  together  for  each 
individual  station.  They  were  quite  noisy  and  were  filtered  to  improve  the  signal.  First,  the  mean 
was  removed  from  each  seismogram  and  high-pass  filtered  with  a  corner  frequency  of  0.5  Hz.  The 
trace  was  cut  to  +/-  100  s  of  the  expected  S  wave  to  detect  any  amplitude  or  frequency  changes 
indicative  of  an  S  arrival.  Only  one  explosion  produced  any  obvious  arrivals,  and  it  was  recorded 
on^only  one  station  (discussed  later).  A  variety  of  different  filters  were  applied  to  the  seismograms 
to  improve  the  possibility  of  detecting  any  S  phases.  A  band-pass  filter  with  corner  frequencies 
at  0.01  Hz  and  0.1  Hz  gave  the  best  results.  By  removing  any  frequencies  with  periods  greater 
than  100  s,  the  long  period  effects  of  the  seismometer  were  eliminated.  These  were  plotted  using 
SAC  in  the  same  manner  as  the  high-pass  filter  seismograms.  Next,  the  components  were  aligned 
according  to  the  expected  S  arrival  for  each  explosion.  The  GMT  (General  Mapping  Tool)  was 
used  to  plot  the  components  versus  epicentral  distance.  The  1988  event  was  recorded  at  only  one 
station.  From  the  GMT  plots,  a  frequency  and  amplitude  variation  can  be  seen  for  a  few  stations 
on  most  components  of  the  8-16-90,  5-21-92,  5-15-95,  and  6-8-96  events  (Figs.  6.1-6.12)  at  the 
expected  S  wave  arrival  at  100  s. 

Discussion 

The  1990  event  had  only  one  visible  arrival  at  the  expected  time-on  the  east  component  of 
station  ARU  (Arti,  Russia).  The  1992  event  had  several  possible  S  waves.  They  were  seen  on  all 
components  at  stations  MAJO  (Matsushiro,  Japan)  and  OBN  (Obninsk,  Russia),  on  the  north  and 
vertical  components  at  KIV  (Kislovodsk,  Russia),  on  the  north  and  east  components  at  KONO 
(Kongsberg,  Norway),  and  on  the  east  component  of  COL  (College  Outpost,  Alaska).  Possible  S 
wave  arrivals  are  seen  on  all  components  of  the  1995  explosion  at  stations  ABKT  (Alibek,  Turk¬ 
menistan),  CHTO  (Chiang  Mai,  China),  and  ARU.  There  also  appears  to  be  a  disturbance  at 
KONO.  However,  it  is  doubtful  it  is  an  S  wave  due  to  the  duration  of  the  signal.  There  also  ap¬ 
pears  to  be  a  phase  arrival  at  TLY  (Talaya,  Russia).  This  is  most  likely  due  to  source  effects.  This 
station  is  out  of  the  epicentral  limit  set  for  S  wave  arrivals  and  was  only  included  because  of  its 
close  proximity  to  Lop  Nor.  Lastly,  the  most  recent  event  produced  possible  S  waves  at  ANTO 
(Ankara,  Turkey),  CHTO  and  ARU  on  the  north  component. 

Analysis 

Because  the  1992  event  is  the  largest  in  magnitude  (M5  =  6.5),  it  was  chosen  as  the  explosion 
to  analyse  further.  At  MAJO,  the  band-pass  filter  produces  a  large  arrival  on  all  three  components 
(Fig.  6.13).  However,  the  PcS  phase  arrives  only  7.3  s  after  the  S  wave,  and  the  filter  may  have 
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combined  the  two  phases.  The  original  seismogram  was  studied  to  determine  if  the  phases  can  be 
seen  separately.  They  are  not  discernible  in  the  original  data,  but  they  are  very  apparent  on  the 
high-pass  filtered  signal  with  a  corner  frequency  of  2  Hz  and  rotated  into  radial  and  transverse 
components  (Fig.  6.14).  There  are  several  puzzling  aspects  to  this  record.  Firstly,  if  the  phases  are 
indeed  S  and  PcS,  then  the  S  wave  is  about  15  seconds  early.  Secondly,  at  an  epicentral  distance  of 
aliiiost  40  degrees,  there  should  be  very  little  vertical  S  phase  energy.  Lastly,  the  PcS  phase  should 
be  almost  completely  radial.  The  substantial  amount  of  energy  on  the  tangential  component  (even 
more  than  on  the  radial  component!)  at  the  expected  arrival  time  of  PcS  requires  explanation.  The 
arrivals  are  not  smface  wave  Love  or  Raleigh  phases  since  these  waves  arrive  around  100s  and  200s 
after  the  expected  S  wave,  respectively. 

The  records  at  OBN  also  show  significant  energy  on  all  three  components  (Fig.  6.15).  The  band¬ 
pass  filter  was  narrowed  to  0.04  Hz  and  0.1  Hz  in  order  to  remove  long  period  effects  completely. 
Again,  the  PcS  phase  arrives  fairly  soon  (30  s)  after  the  S  wave,  and  the  filter  may  have  combined 
the  two  phases.  The  original  seismogram  on  rotated  co-ordinates  shows  the  two  separate  phases 
(Fig.  16).  Again  there  is  substantial  energy  on  the  tangential  component  for  the  PcS  wave. 
In  order  to  determine  if  the  transverse  energy  is  the  result  of  anisotropy,  the  seismogram  was 
cut  into  two  plots  5  s  long  containing  the  S  wave  and  the  PcS  wave  separately,  and  the  north- 
south/east-west  particle  motion  was  plotted  (Figs.  6.17-6.18).  For  waves  passing  through  a  single 
layer  azimuthally  anisotropic  medium,  these  plots  would  be  elliptical.  Fig.  18  does  show  some 
elliptical  motion,  but  Fig.  17  does  not.  A  possible  explanation  for  this  is  that  the  PcS  wave 
passes  through  an  anisotropic  medium  beneath  the  station.  The  S  wave,  on  the  other  hand,  may 
pass  through  anisotropic  material  beneath  the  test  site  and  the  station  resulting  in  two  or  more 
layer  anisotropy  effects  which  as  we  showed  in  the  previous  section  can  give  complicated  tangential 
motion.  Alternatively  the  tangential  motion  could  be  a  combination  of  anisotropy  and  scattering 
The  possible  anisotropy  and  its  parameters  beneath  OBN  are  currently  being  analysed. 

The  anisotropy  study  that  is  being  done  for  OBN  will  be  repeated  for  MAJO.  If  there  are 
anisotropy  effects,  it  could  explain  the  PcS  tangential  motion.  Further  analysis  of  PcS  waves  at 
KIV,  COL,  and  KONO  needs  to  be  done  before  it  can  be  determined  whether  the  study  can  be 
undertaken  at  these  stations  also. 

It  is  still  too  early  to  determine  the  origins  of  the  arrivals  in  the  data  set  and  more  work  is 
required  to  confirm  that  they  are  indeed  S  waves.  The  vertical  energy  must  be  explained.  In 
addition,  the  transverse  energy  seen  at  many  stations  for  the  four  explosions  needs  to  be  explained. 
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The  anisotropy  study  can  be  repeated  for  any  station  that  shows  PcS  waves.  With  these  PcS 
phases,  the  anisotropy  beneath  the  station  can  be  determined  easily  and  compared  with  that  of 
the  previous  section  which  relies  on  using  SKS  phases  from  multiple  earthquakes.  Then  it  would 
be  possible  to  recover  the  anisotropy  beneath  the  test  site.  Such  steps  are  necessary  in  order  to 
recover  the  physical  mechanism  giving  rise  to  teleseismic  S  from  nuclear  explosions. 

•s, 

7.  CONCLUSTONS  AND  RELEVANCE  TO  MONITORING 
THE  COMPREHENSIVE  TEST  BAN  TREATY 

In  sections  1-4  we  described  the  Baikal  rift  zone  experiment  and  the  use  of  teleseismic  body  waves 
to  determine  that  the  asthenosphere  uplifts  in  broad  region  hundreds  of  km  wide.  The  upwarp 
is  associated  with  P  wave  and  S  waves  delays  the  latter  being  three  times  the  former.  It  is  also 
associated  with  strong  attenuation  of  high  frequencies.  We  suspect  this  mantle  anomaly  has  resulted 
from  small  scale  convective  motion  in  which  hot  material  has  upwelled  beneath  the  rift.  Pervasive 
anisotropy  is  observed  a  phenomenon  expected  to  accompany  finite  shearing  of  convecting  mantle 
material.  In  section  5  we  compared  normal  mode  synthetics  with  SKS  waves  and  find  that  nearly 
every  station  of  the  GSN  is  underlain  by  anisotropic  mantle.  Detailed  study  of  ANMO  and  PAS 
reveals  that  about  one  third  of  the  arrivals  have  more  complicated  waveforms  than  expected  from 
a  single  anisotropic  layer  and  are  better  explained  by  a  second  layer  of  anisotropy.  However  the 
second  layer  lies  remote  from  the  station  perhaps  at  the  core-mantle  boundary  and  varies  from 
event  to  event.  Alternatively  these  waveform  complications  could  be  due  to  scattering. 

In  section  6  we  study  S  waves  generated  from  Nuclear  explosions  at  Lop  Nor.  Nuclear  explosions 
are  not  normally  associated  with  S  waves.  However  Lop  Nor  nuclear  explosions  generate  observable 
arrivals  at  the  times  of  teleseismic  S  that  could  have  been  generated  at  the  source  by  cavity  collapse, 
tectonic  stress  release,  or  P-S  conversions.  P-S  conversions  at  the  core-mantle  interface,  PcS,  can 
explain  some  of  the  arrivals.  Elliptical  particle  motion  of  these  arrivals  is  explicable  as  a  result  of 
receiver-side  anisotropy.  S-waves  generated  at  the  source  would  be  affected  both  by  receiver-side 
anisotropy  as  well  as  source-side  anisotropy.  Recognition  of  such  effects  is  needed  to  reconstruct  the 
original  S  wave  polarisation  and  thereby  compare  the  various  candidate  mechanisms  of  generation 
of  nuclear  explosion  S.  The  ultimate  objective  is  to  ascertain  if  the  particle  motion  is  aligned  in 
the  direction  of  inferred  tectonic  stress  release,  or  corresponds  to  P-S  conversions  that  would  be 
radially  directed. 
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Most  teleseismic  location  schemes  assume  the  Earth  is  laterally  homogeneous.  However  in  the 
Baikal  rift  zone  a  localised  delay  of  1  second  occurs  flanked  by  early  arrivals  of  0.3  seconds  either 
side.  These  delays  are  dependent  on  azimuth  and  so  will  effect  event  location  in  this  area.  Similar 
delays  and  structures  have  been  observed  in  the  Rio  Grande  Rift,  and  the  Basin  and  Range  rift 
province  of  North  America  and  on  the  plateau  uplift  surrounding  the  east  African  rift.  All  these 
regions  correspond  to  negative  Bouguer  anomalies  and  near-normal  crustal  thicknesses,  except  in 
the  rift  zones  themselves  where  crustal  thinning  of  about  10  km  is  observed.  Thus  the  delaying 
structures  are  probably  in  the  mantle  and  this  is  corroborated  by  parallax  shifts  in  the  patterns 
of  residuals  as  a  function  of  incident  direction.  We  presume  such  mantle  anomalies  exist  beneath 
many  other  continental  regions  throughout  the  world  and  this  work  suggests  a  strategy  to  use 
gravity  in  conjunction  with  known  crustal  thickness  for  their  identification. 

Probably  the  greatest  implication  of  low  Q  low  velocity  mantle  under  continental  tectonic  zones 
are  the  effects  of  attenuation.  Ringdal  et  al.,  (1992)  find  that  the  yield  Y  in  k.t.  for  nuclear 
explosions  at  the  Shagan  River  test  site,  Russia,  is  related  to  body  wave  magnitude  by 

mb  =  0.75  log  Y  +  4.43 

for  body  waves  in  the  2-4  Hz  frequency  range.  We  can  estimate  the  underestimate  of  yield  for 
regions  of  such  as  a  rift  zone  owing  to  upper  mantle  attenuation.  Let  the  contrast  in  t*  be  6t*  then 
the  ratio  of  yields  in  the  two  zones  is  given  by. 


YifY2  = 

For  example  a  6t*  of  0.1  results  in  a  reduction  of  estimated  yield  of  0.28  of  the  value  had  the 
explosion  been  detonated  outside  of  the  rift  assuming  body  waves  of  3  Hz  are  used  in  the  estimate. 
This  places  extra  requirements  on  techniques  to  monitor  suspected  nuclear  tests  in  continental  rift 
zones  or  regions  above  upwelling  hot  mantle. 

Convection  in  the  mantle  beneath  continents  brings  hot  material  from  depth  which  has  the  effect 
of  lowering  P  and  S  wave  velocities  as  well  as  Q,  and  inducing  anisotropy.  As  a  result,  estimates 
of  yield  based  on  body  waves  could  be  as  much  as  a  factor  of  three  too  low  for  events  detonated 
above  such  an  upwelling  current.  Global  tomography  maps  (Woodhouse,  1995  in  Press,  personal 
communication)  indicate  such  low  velocity  regions  are  common  beneath  orogenic  and  rift  zones  on 
the  continents.  Accurate  location  of  nuclear  events  will  have  to  take  both  lateral  heterogeneity  in 
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velocity  from  such  thermal  effects  as  well  as  that  from  anisotropy  induced  by  the  flow.  Seismic 
anisotropy  can  result  in  converted  phases,  and  has  strong  effect  on  surface  waves. 

Our  work  indicates  that  modelling  of  seismic  waves  using  standard  Earth  models  such  as  PREM 
requires  the  following  corrections  to  fit  the  data. 

1.  Inclusion  of  localised  S  and  P  wave  travel  time  anomalies. 

-'2.  Recognition  of  localised  attenuation  of  high  frequency  P  and  S. 

3.  Removal  of  single  layer  sphtting  of  S  waves. 

4.  Corrections  for  deep  mantle  second  layer  splitting  or  scattering  from  lateral  heterogeneity. 
This  work  provides  a  strategy  for  extension  to  other  orogenic  tectonically  active  continental 

regions  where  global  models  have  insufficient  resolution,  but  geophysical  signatures  such  as  gravity, 
topography,  heat  flow  as  well  as  tectonic  history  could  be  used  to  infer  lateral  variations  in  the 
upper  mantle,  thereby  extending  the  global  models.  The  resulting  model  of  the  continents  would 
then  serve  as  a  basis  for  inverting  seismograms  from  potential  test  sites  for  location  and  yield,  when 
monitoring  a  CTBT. 
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9.  FIGURE  CAPTIONS 


Figure  2.1:  USA-USSR  Baikal  1991  Seismic  Array  Study  UCLA  station  Locations. 

Figure  2.2:  Geographic  location  of  earthquakes  with  mb  >  5.0  occurred  during  July  7  and  Oct 
5,  1991  (Source:  NOAA). 

Figure  2.3:  Comparison  of  teleseisms  and  spectra  from  a  Northern  California  event  from  the  4 
broadband  stations. 

Figure  2.4:  A  seismic  section  from  a  teleseism  from  1  Hz  vertical  sensors,  on  data  stream  1,  i.e. 
the  continuous  stream. 

Figure  3.1:  USA-USSR  Baikal  1992  seismic  array  study  UCLA  station  locations. 

Figure  3.2:  Geographic  location  of  earthquakes  with  mb  >  5.0  recorded  by  UCLA  stations 
during  June  11  and  Oct  5,  1992.  (Source:  NOAA) 

Figure  3.3:  Seismograms  for  teleseismic  events  recorded  at  profile  stations. 

Figure  4.1a:  Polynomial  surface  determined  by  downward  projection  of  the  residual  patterns 
of  Figure  3.  The  X  axis  is  taken  along  the  line  connecting  the  two  stations  at  the  end 
of  the  profile,  and  the  Y  axis  is  perpendicular  to  X  and  pointing  to  the  Northeast.  The 
origin  of  the  co-ordinate  system  is  at  the  middle  point  of  the  lake,  along  the  profile.  The 
contours  on  the  top  and  bottom  planes  are  depth  of  the  asthenosphere  in  kilometres.  The 
shaded  area  on  the  top  plane  represents  Lake  Baikal.  The  uncertainty  of  the  depth  of  the 
polynomial  surface  is  ±17  km. 
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Figure  4.1b:  Regional  Bouguer  gravity  anomaly  and  cross  section  of  Figure  4.1  along  the  profile. 
Crustal  thicknesses  are  taken  from  deep  seismic  sounding.  The  averaged  misfit  between 
the  calculated  gravity  anomaly  (solid  line  in  Figure  a)  from  our  deep  structure  (Figure  b) 
and  measured  regional  anomaly  (circles  in  Figure  a)  is  about  1.0  mGal.  The  best  fitting 
asthenosphere/lithosphere  density  contrast  is  -0.6%. 

^Figure  4.2a:  Mean  travel  time  residual  curves  from  the  1992  array  event  groups.  Eax;h  vertical 
unit  represents  one  second.  Each  horizontal  grid  line  is  the  zero  line  of  a  group,  with 
group  name  written  on  the  right  side  of  the  diagram.  The  two  vertical  curves  connect  the 
location  of  the  two  minima  on  the  residual  curves.  The  insert  on  the  right  shows  average 
arrival  direction  of  each  event  group. 

Figure  4.2b:  Isotropic  velocity  model  determined  from  Bayesian  non-linear  inversion  of  the 
travel  time  residuals  using  a  polynomial  to  represent  the  undulations  of  the  litho- 
sphere/asthenosphere  boundary  that  best  explains  the  travel  time  variation. 

Figure  4.2c:  Mean  S-wave  travel  time  residuals  from  4  events  from  the  South  relative  to  the 
IASP91  earth  model. 

Figure  4.3a:  Map  showing  the  region  with  the  stations  of  the  seismic  array  extending  from  the 
Siberian  platform  across  the  Baikal  rift  zone  into  Mongolia  and  results  of  SKS  splitting 
measurements..  Stations  with  well-defined  measurements  are  represented  by  single  circles 
with  size  proportional  to  the  splitting.  The  line  drawn  through  each  circle  gives  the  fast 
polarisation  direction.  Those  with  2  inconsistent  results  are  plotted  as  double  circles. 
The  shorter  and  longer  line  indicates  the  fast  direction  of  the  smaller  and  larger  circle, 
respectively.  Stations  represented  by  squares  are  those  on  which  anisotropy  effects  cannot 
be  clearly  observed. 

Figure  4.3b:  Normalised  rose  diagrams  of  SKS  fast  polarisation  directions  for  each  of  the  6  areas 
in  the  BRZ  and  adjacent  areas.  Null  measurements  were  excluded  in  the  rose  diagrams. 
The  sector  width  of  the  rose  diagrams  is  20°. 

Figure  4.4:  Results  of  St*  measurements.  The  dots  are  averaged  values  over  all  the  events  for 
St*  (Figure  A)  and  St^  (Figure  B),  and  the  solid  lines  are  the  results  of  spatial  averaging 
in  windows  of  200  km  wide,  moving  at  steps  of  10  km.  The  thickness  of  the  shaded  area 
is  2  times  of  the  standard  error  of  the  mean. 

Figure  4.5:  Relations  between  St*  and  St*.  Note  that  the  values  used  are  from  spatial  averaging 
of  the  mean  values  (the  solid  curves  in  Figure  4.5).  The  best  fitting  results  and  the  cross 
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correlation  coefficient  are  shown  in  the  figure. 

Figure  4.6:  A  Lambert  projection  map  showing  locations  of  events  (numbered  dots)  used  in 
lower  mantle  structure  study,  stations  (open  triangles),  and  ScS  or  PcP  bounce  points 
(filled  triangles). 

Figure  4.7:  Results  of  travel  time  residual  measurements  for  Event  1  plotted  against  latitudes 
of  ScS  bounce  points. 

Figure  4.8:  PcP-P,  ScP-P,  and  ScS-S  residuals  for  events  1,  2,  and  3  plotted  at  latitudes  of  ScS 
bounce  points. 

Figure  5.1:  SKS  waveforms  (dotted  lines)  for  GSN  station  ANMO  compared  to  normal  mode 
synthetics  (smooth  curves,  Woodhouse,  1988).  Radial  seismograms  are  on  the  left.  Az¬ 
imuthal  components  modelled  as  a  double  layer  of  anisotropy  are  in  the  centre  and  as 
a  single  layer  on  the  right.  On  the  far  right  numbers  give;  event  number,  azimuth,  fast 
directions  (anticlockwise  from  north),  and  delay  in  seconds  for  the  upper  and  lower  layers, 
and  CMT  event  name.  Note  events  21,23,24,29  and  30  are  adequately  described  by  a 
single  layer  of  splitting.  However  event  27  and  possibly  28  are  better  described  by  the 
double  layer. 

Figure  5.2:  SKS  waveforms  from  the  64  analysed  at  ANMO  for  which  a  2  layer  anisotropy  is  a 
good  or  better  fit  than  one  layer.  The  fact  that  the  second  layer  parameters  are  different 
and  incident  rays  close  to  these  from  other  events  did  not  exhibit  second  layer  effects  leads 
us  to  believe  these  effects  are  coming  from  deep  in  the  mantle;  quite  possibly  in  the  D” 
region  of  the  lowermost  mantle. 

Figure  6.1  -  6.12:  East,  north  and  vertical  components  of  Lop  Nuclear  Explosions  plotted  as  a 
function  of  epicentral  distance,  and  centred  on  the  expected  time  of  teleseismic  S  were  it 
generated  by  the  explosions. 

Figure  6.13:  Possible  S  wave  arrivals  at  MAJO.  Interference  between  the  S  and  PcS  waves  gives 
the  apparent  long  period  structure  of  the  filtered  seismograms. 

Figure  6.14:  High  pass  filtered  version  of  the  record  shown  in  Figure  6.13,  showing  candidate  S 
and  PcS  arrivals. 

Figure  6.15:  Possible  S  wave  arrivals  at  OBN.  Interference  between  the  S  and  PcS  waves  gives 
the  apparent  long  period  structure  of  the  filtered  seismograms. 

Figure  6.16:  High  pass  filtered  version  of  the  record  shown  in  Figure  6.15,  showing  candidate  S 
and  PcS  arrivals 
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Figure  6.17:  Particle  motion  of  the  S  motion  at  OBN  showing  complicated  polarisation  effects 
that  may  be  due  to  a  combination  of  source  side  and  receiver  side  anisotropy  and/or 
scattering. 

Figure  6.18:  Particle  motion  of  the  PcS  wave  at  OBN.  The  simple  shape  is  consistent  with 
receiver-side  scattering. 
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Figure  5.1:  SKS  waveforms  (dotted  lines)  for  GSN  station  ANMO  compared  to  normal  mode 
synthetics  (smooth  curves,  Woodhouse,  1988).  Radial  seismograms  are  on  the  left.  Azimuthal 
components  modelled  as  a  double  layer  of  anisotropy  are  in  the  centre  and  as  a  single  layer  on  the 
right.  On  the  far  right  numbers  give;  event  number,  azimuth,  fast  directions  (anticlockwise  from 
north,  and  delay  in  seconds  for  the  upper  and  lower  layers,  and  CMT  event  name.  Note  events 
21,23,24,29  and  30  are  adequately  described  by  a  single  layer  of  splitting.  However  event  27  and 
possibly  28  are  better  described  by  the  double  layer. 
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Figure  5.2:  SKS  waveforms  from  the  64  analysed  at  ANMO  for  which  a  2  layer  anisotropy  is  a 
good  or  better  fit  than  one  layer.  The  fact  that  the  second  layer  parameters  are  different  and 
incident  rays  close  to  these  did  not  exhibit  second  layer  effects  leads  us  to  believe  these  effects  are 
coming  from  deep  in  the  mantle;  quite  possibly  in  the  D”  region  of  the  lowermost  mantle. 
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Figure  1.  East  component  of  seismograms  for  explosion  8/16/90  plotted  according  to 
epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  A  disturbance  is  noticeable  at  station  ARU. 
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Figure  2.  North  component  of  seismograms  for  explosion  8/16/90  plotted  according  to 
epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicator  ihe  expected  S  wave  arrival.  No  disturbances  are  noticeable. 
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Figure  3.  Vertical  component  of  seismograms  for  explosion  8/16/90  plotted  according 
to  epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  No  disturbances  are  noticeable. 
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Figure  4.  East  component  of  seismograms  for  explosion  5/21  /92  plotted  according  to 
epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  A  disturbance  is  noticeable  at  stations  COL 
KONO,  MAJO,  and  OBN. 
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Figure  5.  North  component  of  seismograms  for  explosion  5/21/92  plotted  according  to 
epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  A  disturbance  is  noticeable  at  stations  KONO, 
MAJO,  and  OBN,  and  KIV. 
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Figure  6.  Vertical  component  of  seismograms  for  explosion  5/21/92  plotted  according 
to  epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  A  disturbance  is  noticeable  at  stations  MAJO, 
and  OBN,  and  KIV. 
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Figure  7.  East  component  of  seismograms  for  explosion  5/15/95  plotted  according  to 
epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  bandpass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  A  disturbance  is  noticeable  at  stations  KONO, 
ABKT,  CHTO,  ARU,  and  TLY. 
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Figure  8.  North  component  of  seismograms  for  explosion  5/15/95  plotted  according  to 
epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Etz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  A  disturbance  is  noticeable  at  stations  KONO, 
ABKT,  CHTO,  ARU,  and  TLY. 
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Figure  9.  Vertical  component  of  seismograms  for  explosion  5/15/95  plotted  according 
to  epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  A  disturbance  is  noticeable  at  stations  KONO 
ABKT,  CHTO,  ARU,  and  TLY. 
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Figure  10.  East  component  of  seismograms  for  explosion  6/8/96  plotted  according  to 
epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicat  's  the  expected  S  wave  arrival.  No  disturbances  are  noticeable. 
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Figure  11.  North  component  of  seismograms  for  explosion  6/8/96  plotted  according  to 
epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
filtered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  A  disturbances  is  noticeable  at  stations  ANTO, 
CHTO,  and  ARU. 
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Figure  12.  Vertical  component  of  seismograms  for  explosion  6/8/96  plotted  according 
to  epicentral  distance.  The  seismograms  had  the  mean  removed  and  were  band-pass 
hltered  with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.  The  line  down  the  middle 
indicates  the  expected  S  wave  arrival.  No  disturbances  are  noticeable. 
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Figure  6.13 


Figure  13.  East,  north,  and  vertical  components  of  explosion  5/21/92  recorded  at 
MAJO.  The  mean  was  removed  from  the  seismogram  and  it  was  band-pass  filtered 
with  corner  frequencies  of  0.01  Hz  and  0.1  Hz.,  This  filter  may  have  combined  the  S 
and  PcS  phases.  The  vertical  line  indicates  the  expected  S  wave  arrival.  The  bottom 
axis  is  the  time  after  the  arrival  in  seconds. 
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appears  to  be  none  for  this  phase. 
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Figure  6.17 


OBN  BHE 


^xp^ion  “I™  °‘  e^P<=cted  PcS  wave  of 

The  elliptical  shape'So'^  *th^*  would  indicate  anisotropy. 

P  cal  shape  shows  there  may  be  some  for  this  phase. 


Figure  6.18. 


Table  2.1.  Baikal  1991  Station  Information 


Station 

Station 

Coordinates 

Ele. 

Duration 

Sensor 

No. 

Name 

Lat 

Lon 

(m) 

Type 

01 

Novoselova 

53.620°N 

102.645°E 

550 

07/12-10/02 

L4C+Guralp 

02 

Suhoy-Zaglik 

53.246°N 

103.767°E 

468 

07/11-10/03 

L4C 

03 

Kukunut 

52.984°N 

104.714°E 

610 

07/09-10/03 

L4C 

04 

Turgenevka 

53.034°N 

105.648°E 

820 

07/08-10/05 

L4C 

05 

Tyrgan 

52.767°N 

106.345°E 

910 

07/07-10/05 

L4C+Guralp 

06 

Sakhoya 

52.464°N 

107.379°E 

871 

07/19-10/02 

L4C+STS2 

07 

Solonechanyu 

52.264°N 

108.273°E 

830 

07/21-10/02 

L4C+Guralp 

08 

Hasurta 

52.207°N 

109.078°E 

810 

07/22-10/02 

L4C 

09 

Hovinsk 

52.118°N 

110.024'’E 

772 

07/30-10/02 

L4C 

10 

Uldurga 

52.172°N 

110.785°E 

994 

07/30-10/02 

L4C+Guralp 

Table  2.2.  UCLA  Baikal  1991  DAS  Movement  Chart 


Station 

No. 

Start 

Duration 

End 

DAS 

No. 

01 

193:06 

275:10 

369 

02 

192:05 

276:03 

231 

03 

189:03 

276:10 

227 

04 

189:05 

247:08 

241 

252:10 

267:05 

084 

267:06 

278:03 

115 

05 

190:00 

278:01 

338 

06 

200:08 

257:02 

346 

07 

202:00 

263:00 

380 

263:08 

276:00 

241 

08 

203:08 

275:08 

388 

09 

210:01 

275:04 

378 

10 

211:12 

230:22 

084 

231:10 

264:02 

115 

264:03 

275:05 

380 

Table  2.3.  Baikal  1991  UCLA  Reftek  Station  Recording  Parameters 


Data 

Stream  Number 

Sample 

Rate  (sps) 

Recording 
length  (s) 

Data 

Form  (bits) 

Stations 

01 

10 

1800 

16 

ALL 

02 

50 

180 

16 

ALL 

03 

1 

3600 

32 

1,5,7,10 

Table  2.4.  Baikal  1992  Event  Statistics 


Group 

No. 

Mag 

(mb) 

Delta 

(deg) 

Depth 

(km) 

Number  of 
(Events) 

Comments 

1 

>5.0 

00-180 

000-700 

314 

Medium  &  Strong  events 

2 

>5.5 

00-180 

000-700 

54 

Strong  events 

3 

>5.0 

30-085 

000-700 

164 

Tele-events  with  p  as  first  arrivals 

4 

>4.0 

00-030 

000-700 

111 

Significant  regional  &  local  events 

5 

>5.0 

00-180 

100-700 

82 

Medium  &  strong  deep  events 

Table  3.1.  Baikal  1992  Station  Information 


Station 

No. 

Station 

Name 

Coordinates 

Lat  Lon 

Ele. 

(m) 

Duration 

Sensor 

Type 

00 

Bratsk 

55.965°N 

101.410°E 

376 

08/17-10/03 

L4C+STS2 

01 

Pokosnoe 

55.678‘’N 

100.990'’E 

453 

06/27-08/16 

L4C+STS2 

02 

Naratay 

55.052°N 

101.850°E 

346 

06/25-07/21 

S13 

03 

Ust~kada 

54.5 16°N 

102.070°E 

433 

06/23-10/01 

L4C 

04 

Sborniy 

54.193°N 

102.649°E 

624 

06/22-10/01 

S13 

05 

Konovalova 

53.929°N 

102.934°E 

386 

06/21-09/30 

L4C 

06 

Melchituy 

53.649°N 

103.255'’E 

355 

06/20-09/30 

S13+Guralp 

07 

Suhoy-zaglik 

53.243°N 

103.767°E 

719 

06/18-09/29 

L4C 

08 

Zacharovska 

52.993'>N 

103.927“E 

480 

07/02-08/11 

HSIO 

09 

Baroy 

52.778°N 

104.105°E 

518 

06/17-07/17 

L4C 

10 

Goryashina 

52.622°N 

104.234°E 

524 

06/14-09/29 

HSIO 

11 

Patrony 

52.169°N 

104.469°E 

619 

06/13-09/16 

L4C+STS2 

12 

Listvyanka 

51.847°N 

104.893“E 

600 

06/11-10/05 

S13 

13 

Naratay 

55.022°N 

102.055°E 

388 

07/21-10/01 

S13 

14 

Stepanovka 

52.854°N 

103.966°E 

619 

07/17-09/29 

L4C 

15 

Klyuchi-bulak 

55.560°N 

101.803°E 

400 

08/18-10/03 

L22+Guralp 

16 

Patrony2 

52.162°N 

104.464°E 

419 

09/19-10/05 

L4C+STS2 

80 

Suha-bator2 

50.193°N 

106.254°E 

608 

07/17-10/02 

L4C+Guralp 

81 

Suha-bator 

50.242°N 

106.240°E 

750 

06/15-06/17 

L22+STS2 

82 

Yoroo 

49.738°N 

106.202°E 

682 

06/16-07/17 

S13 

83 

Sharingol 

49.288°N 

106.412°E 

974 

06/17-10/02 

L4C 

84 

Hara 

48.931°N 

106.682°E 

949 

06/17-10/02 

S13+Guralp 

85 

Tunhe 

48.383°N 

106.783°E 

1122 

06/21-10/01 

L4C 

86 

Ulan-bator 

47.921°N 

106.954°E 

1284 

06/20-07/21 

S13 

87 

Bayan-suma 

47.209°N 

107.422'’E 

1320 

07/09-09/26 

L22+Guralp 

88 

Sumber-suma 

46.635°N 

107.758°E 

1265 

07/10-09/26 

S13 

89 

Bayan-argalan 

46.115''N 

107.619°E 

1477 

07/10-09/25 

L4C 

90 

Endershil 

45.262°N 

108.260°E 

1232 

07/11-09/25 

L22+STS2 

91 

Yoroo2 

49.747'’N 

106.188°E 

629 

07/17-10/02 

S13 

92 

Ulan-bator2 

47.921°N 

106.954°E 

1284 

07/21-09/28 

S13 

Table  3.2.  UCLA  Baikal  1992  DAS  Movement  Chart 


Station 

No. 

Duration 

Start 

End 

DAS 

No. 

00 

230 

276 

119 

01 

179 

229 

119 

02 

177 

203 

373 

03 

175 

275 

376 

04 

174 

275 

395 

05 

173 

274 

099 

06 

172 

274 

369 

07 

170 

273 

378 

08 

184 

200 

345 

09 

169 

199 

363 

10 

166 

273 

347 

11 

165 

260 

359 

12 

163 

260 

367 

260 

261 

153 

261 

266 

359 

13 

203 

275 

373 

14 

199 

273 

363 

15 

231 

276 

342 

16 

263 

273 

367 

80 

199 

276 

147 

81 

167 

168 

239 

82 

168 

199 

152 

83 

169 

276 

361 

84 

169 

276 

339 

85 

173 

275 

388 

86 

172 

203 

346 

87 

191 

270 

374 

88 

192 

270 

231 

89 

192 

269 

154 

90 

193 

269 

341 

91 

199 

276 

152 

92 

203 

272 

346 

Table  3.3.  Baikal  1992  UCLA  Reftek  Station  Recording  Parameters 


Data 

Stream  Number 

Sample 
Rate  (sps) 

Recording 
length  (s) 

Data 

Form  (bits) 

Stations 

01 

10 

1800 

16 

ALL 

02 

25 

180 

16 

STS2 

50 

180 

16 

Guralp;  short  period 

03 

1 

1800 

32 

Guralp;  STS2 

Table  3.4.  Baikal  1992  Event  Statistics 


Group  Mag  Delta  Depth  Number  of  Comments 

No.  (mb)  (deg)  (km)  (Events) 


1 

>5.0 

000-180 

000-700 

464 

2 

>5.5 

000-180 

000-700 

113 

3 

>5.0 

300-085 

000-700 

217 

4 

>4.0 

000-030 

000-700 

191 

5 

>5.0 

000-180 

100-700 

81 

Medium  &  Strong  events 
Strong  events 

Tele-events  with  p  as  first  arrivals 
Significant  regional  &  local  events 
Medium  &  strong  deep  events 


